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ABSTRACT 

We use deep optical, infrared and radio observations to explore the symbiosis be- 
tween the nuclear activity and galaxy evolution in the southern compact radio source 
PKS1549-79 [z = 0.1523). The optical imaging observations reveal the presence of 
tidal tail features which provide strong evidence that the host galaxy has undergone 
a major merger in the recent past. The merger hypothesis is further supported by the 
detection of a young stellar population, which, on the basis of spectral synthesis mod- 
elling of our deep VLT optical spectra, was formed 50 - 250 Myr ago and makes up a 
significant fraction of the total stellar mass (1 - 30%). Despite the core-jet structure of 
the radio source, which is consistent with the idea that the jet is pointing close to our 
line of sight, our HI 21cm observations reveal significant HI absorption associated with 
both the core and the jet. Moreover the luminous, quasar-like AGN (My < —23.5) 
is highly extinguished (A„ > 4.9) at optical wavelengths and shows many properties 
in common with narrow line Seyfert 1 galaxies (NLSl), including relatively narrow 
permitted lines {FWHM -1940 km s'^), highly blueshifted [OIII]AA5007,4959 hues 
{/SV — 680 km s^^), and evidence that the putative supermassive black hole is ac- 
creting at a high Eddington ratio (0.3 < Lboi/Ledd < H)- The results suggest that 
accretion at high Eddington ratio does not prevent the formation of powerful rela- 
tivistic jets. 

Together, the observations lend strong support to the predictions of some recent 
numerical simulations of galaxy mergers in which the black hole grows rapidly through 
merger-induced accretion following the coalescence of the nuclei of two merging galax- 
ies, and the major growth phase is largely hidden at optical wavelengths by the natal 
gas and dust. Although the models also predict that AGN-driven outflows will even- 
tually remove the gas from the bulge of the host galaxy, the visible warm outflow in 
PKS1549-79 is not currently capable of doing so. However, much of the outflow may 
be hidden by the material obscuring the quasar and/or tied up in hotter or cooler 
phases of the interstellar medium. 

By combining our estimates of the reddening of the quasar with the HI column 
derived from the 21cm radio observations we have also made the first direct estimate 
of the HI spin temperature in the vicinity of a luminous AGN: Tspin > 3000K. 

Key words: galaxies:starburst - galaxies:individual: PKS1549-79 - galax- 
ies:formation - galaxies:interactions - quasars:general 



1 INTRODUCTION 

Powerful radio galaxies have long been considered as use- 
ful signposts to massive structures in the distant universe. 
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However, there is now increasing speculation that AGN and 
associated jet activity are not merely signposts but, in fact, 
intimately linked to the galaxy evolution process. The main 
evidence for such links is twofold. On the one hand the shape 
of the rcdshift evolution in the number density of quasars 
and radio galaxies shows marked similarities to that of the 
star formation history of the universe (Dunlop & Peacock 
1990, Madau et al. 1996), while on the other close correlar 
tions have been found between the masses of the supermas- 
sive black holes that produce the prodigious nuclear activity 
and the properties of the bulges of the host galaxies (e.g. 
Tremaine et al. 2002). 

The correlations between the rcdshift evolution of the 
radio source population and the star formation history of the 
field galaxy population can be explained in terms of hierar- 
chical galaxy evolution models (e.g. Kauffmann & Haenelt 
2000) . It is well known from studies of nearby merging galax- 
ies that stars arc formed as part of the merger process. Nu- 
merical simulations of mergers of gas rich galaxies also show 
that the tidal torques associated with the mergers concen- 
trate gas in the circumnuclear regions of the galaxies (Barnes 
& Hearnquist 1996, Mihos & Hearnquist 1996); some of this 
gas may be available to fuel circumnuclear starburst, quasar 
and jet activity. Therefore it is likely that AGN activity and 
near-nuclear staxbursts in massive galaxies axe a direct con- 
sequence of hierarchical galaxy evolution; AGN may be a 
symptom of the overall galaxy evolution process. 

In this context, radio galaxies are key objects for under- 
standing the link between star formation and AGN activity 
in individual galaxies. Not only is there clear morphological 
evidence that a significant subset of powerful radio galax- 
ies — ~50% of those with strong emission lines — have 
undergone recent mergers (Heckman et al. 1996, Smith & 
Heckman 1989), but many of the radio galaxies with mor- 
phological evidence for mergers also show signs of recent star 
formation (Tadhunter et al. 1996, Arotxaga et al. 2001, Tad- 
hunter ct al. 2002, 2005, Wills et al. 2002, 2004, Johnston 
et al. 2005). Despite these results, considerable uncertainties 
remain about the relationship between star formation and 
AGN/jet activity in radio galaxies. For example, it is not 
certain at what stage of the merger the AGN activity is trig- 
gered, nor whether the star formation is exactly coeval with 
the AGN activity. Indeed, some recent studies of the young 
stellar populations in radio galaxies suggest that the activ- 
ity may be triggered relatively late in the merger sequence, 
a significant time after the major episode of merger-induced 
star formation (Tadhunter et al. 2005, Emonts et al. 2005, 
Johnston et al. 2005). However, the sample sizes are cur- 
rently small and further detailed studies of the young stellar 
populations in powerful radio galaxies are required to put 
our understanding of the link between star formation and 
AGN activity on a firmer footing. 

While the links between star formation and AGN/jet 
activity provide evidence that AGN arc a bi-product of the 
overall hierarchical galaxy evolution process, it has been sug- 
gested that the correlations between the supermassive black 
holes masses and the properties of the bulges of the host 
galaxies are a consequence of the feedback effects of AGN- 
induced outfiows (Silk & Rces 1998, Fabian 1999, di Matteo 
et al. 2005). To date, most of the evidence for such outflows 
has been provided by observations of narrow and broad ab- 
sorption lines observed in the spectra of AGN at UV and 



X-ray wavelengths (see Crenshaw, Kraemer & George 2003 
for a review) . The outflows associated with these absorption 
features clearly have an impact on the near-nuclear environ- 
ments on a scale <10pc, but the impact they have on the 1 
- 10 kpc scale of galaxy bulges is less clear. Although out- 
flows have also been detected in both ionized and neutral 
gas on a kpc-scale in powerful radio galaxies, uncertainties 
remain about the general importance of the outflows and 
their driving mechanism (quasars, jets or starbursts). 

The subject of this paper PKS1549-79 is a key 
object for studies of the links between galaxy evolution and 
AGN activity, because it shows clear evidence of both the 
merging process that may trigger the activity, and the feed- 
back efi'ect of the AGN activity on the surrounding ISM. 

PKS1549-79 is a compact flat spectrum radio source. 
Previously unpublished VLBI radio images of the source 
show that it has the compact core-jet structure expected for 
such sources, consistent with the idea that its jet is point- 
ing close to the line of sight. Remarkably for such a source, 
low-resolution optical spectra fail to reveal the broad per- 
mitted lines and strong non-stellar continuum characteris- 
tic of quasars and BL Lac objects (Tadhunter et al. 2001). 
Therefore this object does not readily fit in with the simplest 
versions of the unified schemes which propose that lines of 
sight close to the radio axis have relatively unobscured views 
of the AGN (Barthel 1989). 

Despite the absence of strong quEisar continuum and 
broad emission line features, optical spectra of the source do 
reveal high ionization forbidden emission lines which suggest 
that the source contains a powerful quasar-like AGN, even 
if this AGN is not directly visible at optical wavelengths. 
However, the high ionization omission lines are unusually 
broad, and blueshifted by ~650 km/s relative to the low 
ionization lines. In addition, this object has HI 21cm absorp- 
tion detected against its radio core, at the same redshift as 
the low ionization omission line system. Further support for 
the idea that the quasar nucleus is obscured by the circum- 
nuclear material in this source is provided by K-band spec- 
troscopic observations that detect a broad Paa emission line 
and strong non-stellar continuum characteristic of quasars 
(Bellamy ct al. 2003). 

Aside from its unusual emission line properties, 
PKS1549-79 belongs to the subset of ~30 - 50% of powerful 
radio galaxies for which young stellar populations malce a 
substantial contribution to the optical continuum emission 
(Tadhunter et al. 2002). Further evidence for star forma- 
tion in this object is provided by its unusually strong far-IR 
emission, which leads to its classification as an ultra lumi- 
nous infrared galaxy (ULIRG: see Sanders & Mirabel 1996 
for a definition) with Lir = 1.6 x 10^^ L© . This evidence 
for recent star formation is consistent with the idea that 
the host galaxy of this source has been involved in a merger 
which has also triggered a major episode of star formation. 

In this paper we present new deep spectroscopic and 
imaging observations of PKS1549-79 at optical, infrared and 
radio wavelengths that allow us to investigate the triggering 
of its AGN activity, characterise its circumnuclear outflows, 
and deduce the nature of its AGN. In this way we can place 
this key object in the context of evolutionary scenarios for 
the origin and evolution of the activity. 

We assume cosmological parameters of Ho = 75 km s~^ 
and qo = 0.0 throughout the paper. For these parameters 
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1 arcsec corresponds to 2.39 kpc at the redshift of the host 
galaxy of PKS1549-79. 



2 OBSERVATIONS AND REDUCTIONS 
2.1 Optical spectroscopy 

Optical spectroscopic observations were taken in 2002 July 

and 2003 September using the EMMI spectrograph in RILD 
mode on the ESO 3.5-m New Technology Telescope (NTT), 
and the FORS2 spectrograph in LSS mode on the ESO Very 
Large Telescope (VLT). Both sets of data are deeper and 
have a higher spectral resolution than those reported in Tad- 
hunter et al. (2001). Details of the observations are given in 
Table 1. 

On the NTT, grisms #5 and #6 were used with the 
MIT/LL mosaic to obtain spectra with wavelength ranges 
3215-6136A and 4977-7565A in the blue and red respectively. 
These configurations resulted in spectral resolutions of 5.0 
- 6.7A in the blue and 4.2 - 5.6A in the red, depending 
on slit width. Red and blue exposures were interleaved in 
an attempt to ensure uniformity in conditions between the 
exposures. To reduce the effects of differential refraction, all 
exposures were taken with the slit aligned along the paral- 
lactic angle. Spectra were taken along two slit positions - 
PA 25 and PA -5. 

The aim of the NTT spectra was to investigate the 
kinematics and physical conditions of the emission line gas 
through modelling of the bright emission lines. However, 
higher signal-to-noise is required to investigate the stellar 
populations in PKS 1549-79. Hence, we have also obtained 
optical spectra with the F0RS2 spectrograph on the VLT on 
2003 September 25. Grisms 600B and 600RI (holographic) 
were used with the MIT CCD to obtain spectra with wave- 
length ranges 3100-6000A and 5000-8200Arespectively. The 
VLT spectra were taken with a 1.3 arcsec slit to ensure the 
entire nuclear region was sampled, and the slit was aligned 
along the parallactic angle (PA 75). The resulting spectral 
resolutions were 6.5A and 7.4A for the blue and red spec- 
tra respectively. As with the NTT observations, the red and 
blue exposures were interleaved. 

For both the NTT and the VLT observations the data 
reduction followed the standard steps of bias subtraction, 
flat fielding, wavelength calibration, atmospheric extinction 
correction, flux calibration, and registration of blue and red 
spectral images. Based on measurements of night sky lines, 
the uncertainty in the wavelength calibration is estimated 
to be smaller than 0.5A for all datasets; the relative flux 
calibration uncertainty, based on the comparison of obser- 
vations of several spectroscopic standard stars, is better than 
±5%. Prior to the analysis of the data, the spectra were cor- 
rected for a Galactic extinction assuming E{B — V) — 0.679 
(Schlegel, Finkbeiner & Davis 1998) and the Seaton (1979) 
extinction law. 



2.2 Optical imaging 

Deep images of PKS1549-79 were obtained using FORSl 
on the VLT with a Gunn r filter in good seeing conditions 
{FWHM = 0.67 arcseonds). The advantage of this filter is 
that, covering the rest wavelength range ~5200 - 6IOOA, it 



Date Exposure PA Slit width Seeing FWHM 



(s) (circsec) (arcsec) 



NTT/EMMI (Optical Spectroscopy) 

12/07/02 3*1200(R) -5 1.0 0.78 

3*1200(B) -5 1.0 0.78 

13/07/02 2*1200(R) 25 1.5 0.78 

2*1200(B) 25 1.5 0.78 

VLT/FORSl (Gunn r Imaging) 

25/07/03 4*300 - - 0.67 

VLT/FORS2 (Optical Spectroscopy) 

25/09/03 3*600(R) 75 1.3 1.2-1.8 

25/09/03 3*1200(B) 75 1.3 1.2-1.8 

AAT/IRIS2 (K-band Spectroscopy) 

4/09/03 24*120 1.0 1.0 

NTT/SOFI (H and K-band Spectroscopy) 

24/08/04 12*120 1.0 1.2 



Table 1. Log of optical and infrared observations for PKS1549- 
79. In the case of the optical spectroscopic observations the seeing 
was estimated from the DIMM monitor estimates at the time of 
the observations, whereas for the imaging observations the seeing 
was measured directly from star images in the field of the source. 

avoids all strong emission lines in the source and is therefore 
dominated by continuum light. During the observations the 
telescope was jittered to four separate positions on the sky 
to help remove cosmetic defects and average out sensitivity 
variations. 

The imaging data were reduced using a combination of 
the STARLINK FIGARO package and IRAF, by first sub- 
tracting the bias, dividing each separate frame by a jittered 
sky flat fleld, then combining the images with IMCOMBINE 
in IRAF using the coordinate information in the image head- 
ers. 

2.3 Near-IR spectroscopy 

In an attempt to improve on the near-IR spectroscopic data 
of Bellamy et al. (2003), which were affected by a faulty 
grism in IRIS2, new near-IR spectra of PKS 1549-79 were 
taken with the Anglo Australian Telescope (A AT), and the 
NTT. Details of the observations are given in Table 1. 

The first data set was obtained on 2003 September 4 
using the IRIS2 instrument on the AAT. The Sapphire 240 
grism was used to cover the K band. In total, twenty four 
exposures of 120s were obtained for PKS 1549-79, with the 
galaxy 'nodded' between two positions on the slit in an 
ABBA pattern. Four exposures of 66s were taken of the AO 
star HIP 77941 at a similar airmass for calibration purposes. 

The second set was obtained on 2004 August 24 with 
the SOFl instrument on the NTT. The low resolution red 
grism was used with the GRF order-sorting filter to cover 
the H and K bands simultaneously. Cloud was present dur- 
ing the night and the observations were not photometric. 12 
exposures of 120s were obtained; the galaxy was again nod- 
ded along the slit but this time with a random jitter every 2 
exposures. Eight exposures of 30s were taken of the AO star 
HIP 79933 at a similar airmass for calibration purposes. 
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The IRIS2 exposures were co-added using median filter- 
ing to remove cosmic rays; the median 'B' image was then 
subtracted from the median 'A' image to remove the night 
sky lines. The standard star frames were combined in a sim- 
ilar way. The resulting galaxy and star frames wore then 
fiat-fielded using a dark-subtracted flat-field taken with the 
same instrumental setup. The night sky lines from the data 
frames themselves were used to determine a two-dimensionaJ 
wavelength calibration, which was then applied to the fi- 
nal co-added frames of the galaxy and the standard star. 
The resuhing spectra cover the range 20320 - 24900A. The 
spectral resolution, measured from the night sky lines, is 
10.9±0.4 A(FWHM). Fits to the night sky lines show that 
the wavelength calibration is accurate to ±0.16 A. The two- 
dimensional spectra were straightened using the APSUM 
package in IRAF, and one-dimensional spectra were ex- 
tracted with FIGARO. 

Night sky lines were also used to derive a wavelength 
calibration for the NTT data. The resulting two-dimensional 
transformation was then applied to all the individual data 
frames and then suitable groups of exposures (because of 
the random jittering) were combined with FIGARO and a 
series of one-dimensional spectra wore extracted and then 
combined together. There was no appreciable tilting of the 
slit image and so the SOFI spectra did not need to be 
straightened before extraction. The SOFI spectrum cov- 
ers the wavelength range 15057 - 25392A. Atmospheric ab- 
sorption dominates between ~18000 and 19500A and the 
thermal background dominates at wavelengths greater than 
~24000A. The spectral resolution, measured from the arc 
lines, is 35.0±0.4 A(FWHM). Fits to the night sky fines 
show that the wavelength calibration is accurate to ± 0.77 

A. 

The standard star exposures were used to remove tel- 
luric features and also to flux calibrate the IRIS2 spectra. 

The magnitudc-to-flux conversion was performed with ref- 
erence to Bessell, CastcUi & Plcz (1998), while telluric cor- 
rection (and flux calibration for the IRIS2 data) was per- 
formed by assuming the intrinsic SED of the AO stars to be 
that of a black body at T = 9480K. For hot stars this is a 
good approximation, particularly in the near-IR end of the 
black body spectrum. A rmmber of Brackett lines are visi- 
ble in the spectra of AO stars, particularly in the H band, 
but these were masked out before the telluric correction was 
performed. 

The SOFI data suffered unquantifiablc slit-losses due 
to the telescope nodding axis not being aligned precisely 
with the slit. For this reason no attempt was made to make 
an absolute flux calibration of the SOFI spectrum. Instead 
the spectrum was scaled to match the IRIS2 spectrum in 
the K-band and the spectrum of HIP 79933 was used only 
to remove telluric features and obtain the correct spectral 
shape. 

2.4 Radio observations 

Continuum VLBI observations of PKS 1549-79 were made at 
both 2.3 and 8.4GHz with the SHEVE network (Preston et 
al. 1993, Jaunccy ct al. 1994) between November 1988 and 
March 1992, using the MK2 recording system (Clark 1973). 
Full details of the antennae used for SHEVE observations 
are given in Tzioumis et al. (2002). The tapes recorded at 



each station were processed at the Caltoch-JPL Block2 cor- 
relator in Pasadena, and the NRAO AIPS package was used 
for global fringe fitting. The amplitude calibration and edit- 
ing of the visibilities were subsequently performed using the 
Caltech VLBI package. Self-calibration and imaging of the 
data were performed using DIFMAP (Shepherd et al. 1994) 
and the AIPS packages. At 2.3GHz the restoring beams were 
7.3x2.7 mas (PA4.7) and 8.9x7.2 mas (PA9.8) with and 
without the Harteebeesthoek antenna respectively, whereas 
at 8.4 GHz the restoring beam was 6.7x5.2 mas (PA73) 
without the Hartebeesthoek antenna. The r.m.s. noise in 
the final images is 2.7 mjy beam~^ at 2.3 GHz, and 3.8 
mjy beam~^ at 8.4 GHz. 

In order to investigate the spatial distribution of the 
HI 21cm absorption previously detected by Morganti et al. 
(2001), HI VLBI observations were carried out using the 
Long Baseline Array (LBA) at the redshifted HI frequency of 
1645 MHz. At this relatively low frequency we could use only 
three telescopes: Parkes (64m), ATCA (Tied Array, 5x22m) 
and Mopra (22m). We used a 16MHz bandwidth in each 
circular polarization and 1024 channels. The final data cube 
(after Banning smoothing) has a spectral resolution of ~ 9 
km s~^ , and a noise per spectral channel of 5.3 mJy beam~^, 
while the continuum map derived from these data has a noise 
of 16 mJy beam"^. 

The editing part of the calibration of the 21cm line data 
was done in AIPS and then the data were transferred to 
MIRIAD (Sault, Tcubcn & Wright 1995) for the bandpass 
calibration. The calibration of the bandpass was done using 
additional calibrations of the strong calibrators PKS1718- 
649 and PKS1934-638. 

Both the line and the continuum image were made us- 
ing uniform weighting and the restoring beam is 47x32mas 
(PA-77.3). As expected, the resolution is lower than that 
obtained in the previous VLBI continuum observations de- 
scribed above. The total cleaned flux in the continuum image 
is 4.75 Jy — close to that obtained from the previous low 
resolution ATCA observations. 



3 RESULTS 

3.1 Optical imaging 

Our deep VLT image of PKS 1549-79 is shown in Figure 
1. This reveals a disturbed optical continuum morphology, 
with a high surface brightness elongation extending 4 arc- 
seconds (10 kpc) to the north of the rmcleus, a pair of tidal 
tails (or loop feature observed close to edge-on) extending a 
total distance of 10 arcseonds (24 kpc) to the south, and a 
fainter jet-like feature that extends 12 arcseconds (29 kpc) 
to the NNE of the nucleus. The extended structures con- 
tain several knots of emission, particularly to the SSE and 
NNW. Overall, the structures suggest that PKS 1549-79 has 
undergone a recent merger, which induced star formation in 
the extended halo of the galaxy. 

3.2 Optical spectroscopy 

Three slit positions used for the optical spectroscopy obser- 
vations are shown superimposed on the Gunn r image in 
Figure 1. 
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Figure 1. Deep VLT image of PKS1549-79 tciken witli a Gunn r filter that avoids strong emission lines from the galaxy. The slit positions 
used for the long-slit spectroscopic observations are shown for reference. North is to the top and east to the left. 



3.2.1 Optical continuum modelling 

Accurate modelling of the continuum is required not only 
to deduce the properties of any young stellar populations 
present, but also as a prerequisite for accurate measurement 
of the emission lines, following subtraction of the best-fit 
continuum model. In this section we concentrate on the VLT 
spectrum extracted for the nuclear region because this has 
higher S/N in the continuum than either of the NTT spectra. 

Prior to continuum modelling, the level of the nebular 
continuum was deduced following the procedure outlined in 
Dickson et al. (1995) and Holt et al. (2003), based on double 
Gaussian fits to the H/3 emission line feature. This nebular 
continuum was then subtracted from the data. The major 
uncertainty in this procedure is the reddening in the nebular 
contirmum. Unfortunately, estimates of the reddening using 
the Balmer lines are highly uncertain in this object (see sec- 
tion 3.2.3). To consider the effects of the reddening of the 
nebular continuum on the results of the continuum mod- 



elling, we carried out the continuum modelling twice with 
two extreme sets of assumptions about nebular continuum. 

First we assumed zero reddening for the nebular continuum. 
Second we assumed the max;imum reddening allowed by the 
Balmer line ratios. 

The technique used to model the continuum is described 
in detail in Tadhunter et al. (2002, 2005). To summarise, the 
continuum flux is measured in several wavelength bins (in 
this case 39) chosen to avoid emission lines, cosmetic de- 
fects and poorly-corrected atmospheric absorption features. 
A minimum technique is then used to compare the mea- 
sured fluxes with the results of spectral synthesis models. 
Models comprising various combinations of an old stellar 
population (12.5 Gyr), young stellar populations (0.01 - 
2.5 Gyr), reddened young stellar populations (0.01 - 2.5 Gyr, 
0.0 < E{B — V)< 2.0), and a power-law are compared with 
the data. For the spectral synthesis results we used the most 
recent set of high resolution models from Bruzual and Char- 
lot (2003) for an instantaneous burst, a Salpeter IMF and 
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solar metallicity. For the reddened young stellar populations 
(YSP) we used the Seaton (1979) extinction law to redden 
the YSP spectra from the spectral synthesis models. We 
also Eissumed a uniform uncertainty of ±5% for all the mea- 
sured fluxes, consistent with the estimated uncertainty of 
the relative flux calibration. Based on our experience of this 
modelling technique, models with Xred < 1 can be deemed 
to provide a good fit to the data (see Tadhunter et al. 2005 
for discussion). 

In our first attempt to model the nebular-subtracted 
optical continuum, we considered the full wavelength range 
covered by the spectra. It became clear that no combination 
of old and young stellar populations (reddened or otherwise) 
could fit the relatively steep rise in the spectrum to the red 
at wavelengths longer than ~6500A, and that a red power- 
law is required to provide an adequate fit in that spectral 
range. Note that the steep rise to the red is not due to 
any problem with the reduction and calibration of the data, 
since it is also present in our NTT spectra. Figure 2 shows 
the best fitting model that comprises a 12.5 Gyr old stellar 
population, a 0.05 Gyr young stellar population reddened by 
E{B — V) = 0.4, and a power-law. The level of the power- 
law component is consistent with the extrapolation of the 
reddened quasar continuum detected at near-IR wavelengths 
(see section 3.3). Indeed the existence of the red power- law 
provides independent evidence for the partially extinguished 
quasar first detected by Bellamy et al. (2003). 

In order to avoid the added uncertainly that the in- 
clusion of a power-law introduces into the modelling of the 
continuum, and to deduce an accurate age for the YSP, we 
then excluded wavelengths A > 6OOOA, from the model fits. 
In the case of maximum reddening for the subtracted neb- 
ular continuum, it is clear that the data are well fitted by 
models that include a combination of a 12.5 Gyr old stellar 
population and a YSP with age in the range 0.04 - 0.1 Gyr 
and reddening in the range 0.0 < E{B — V) < 0.75, with a 
best fit obtained for a YSP age of 0.05 Gyr and reddening 
of E{B — V) = 0.6. The results appear to be robust in the 
sense that, if we assume that the nebular continuum suffers 
no reddening, the results are similar, except that the range 
of acceptable models extends to older ages for the YSP (0.05 
- 0.25 Gyr) and the best fitting model has an age of 0.1 Gyr 
and a reddening of E(B-V)~0.5. The results of the SED 
modelling for the restricted wavelength range are shown in 
Figure 3 for the zero reddening in the nebular continuum. 

The next step in the analysis was to check the SED 
fitting results by determining how well the models fit the 
profiles of absorption lines that do not suffer from signifi- 
cant emission line contamination, namely the Call K, the 
G-band, and Mgl features. Reassuringly, the models that fit 
the general SED well also provide the best fits to the ab- 
sorption lines. Figure 4 shows an example of a fit to the 
optical spectrum for a model with a YSP of age 0.1 Gyr and 
reddening E{B — V)= 0.5, assuming zero reddening for the 
nebular continuum. 

To summarise: fits to the optical continuum require a 
young stellar population (YSP) with age in the range 0.04 - 
0.25 Gyr and reddening 0.0 < E{B - V) < 0.8 (depending 
on age) in order to provide an adequate fit to the optical 
SED and stellar absorption features. The YSP has a mass 
2 X 10* < Mysp < 4 X 10^ M0 and makes up 1 — 30% of 
the total stellar mass in the slit (depending on the exact age 



and reddening) . Note the age deduced for the YSP is signif- 
icantly smaller than the 1 Gyr determined by Tadhunter et 
al. (2002) using lower resolution data with a narrower wave- 
length range. Much of this difference is likely to be due to 
the fact that Tadhunter et al. (2002) did not correct their 
spectrum for Galactic extinction prior to the modelling. 

3.2.2 Emission line kinematics 

The long-slit spectra of PKS1549-79 provide useful infor- 
mation on the emission line kinematics in both the nuclear 
and off-nuclear emission line regions. The spectra were mea- 
sured by using the DIPSO package to fit Gaussian profiles 
to the lines. The velocity widths derived from the fits were 
quadratically corrected for the instrumental profiles, and all 
line widths and radial velocity shifts were corrected to the 
rest frame of PKS1549-79. 

The NTT spectra, which were obtained under good see- 
ing conditions, provide the best information about the spa- 
tial extent and kinematics of the extended emission line gas. 
Figure 5 presents flux and radial velocity information de- 
rived from the fits to the emission line profiles. The Ha emis- 
sion line is detected across the full extent of the higher sur- 
face brightness continuum structures detcted in the Gunn r 
continuum image shown in Figure 1, extending 10 arcsec- 
onds to the south and 5 arcscconds to the north of the 
rmclcus. The emission line kinematics are relatively quies- 
cent in the extended regions with small velocity amplitudes 
{AV <150 km s"^) and linewidths (FWHM < 400 km s"^) 
— consistent with a gravitational origin for the gas motions 
(Tadhunter et al. 1989, Baum et al. 1990). 

In contrast to the extended regions, the kinematics in 
the nuclear region provide evidence for more extreme gas 
motions, as already noted by Tadhunter et al. (2001). In 
order to analyse the gas motions in the spatially unre- 
solved nuclear regions we have adopted the approach used 
by Holt et al. (2003) for the radio galaxy PKS1345+12. 
This consists of fitting the minimum number of Gaussian 
components required model each line adequately. For this 
analysis we have concentrated on the VLT data which 
have a higher S/N in the nuclear regions. In the case of 
emission line doublets ([OIII]AA5007,4959, [OI]AA6300,6363, 
[NII]AA6548,6584, [SII]AA6717,6731) the separations, and, 
with the exception of the [SII] blend, the ratios of the 
lines were set by atomic physics, and it was assumed that 
both components of the doublets have the same kinematic 
components. We found that for all the low ionization lines 
[OII]A3727, H/3, [OI]AA6300,6363, [NII]AA6548,6584, and 
[SII]AA6717,6731, a double Gaussian fit comprising narrow 
(FWHM < 400 km s"^) and intermediate width blueshifted 
components (FWHM > 1000 km s^^) provided an adequate 
fit to the emission line profiles. For these low ionization lines 
the narrow component is strongest, and the widths and red- 
shifts of the narrow components are consistent within the 
uncertainties for all the narrow lines; they are also consis- 
tent with the redshift of the Ha emission line detected in the 
extended regions (see Figure 5). The intermediate emission 
line component is relatively weak in the low ionization lines, 
and its width and redshift less well determined based on free 
double Gaussian fits. 

Initial attempts to fit the [OIII]AA5007,4959 lines 
showed that it was possible to fit each of line of the dou- 
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Figure 2. Spectral synthesis model fit to the full wavelength range of the optical spectrum of PKS1549-79. The model comprises a 
0.05Gyr YSP reddened by E{B — V) = 0.4, a 12.5Gyr old stellar population, and a power-law. Note that the power-law component is 
required to fit the steep rise in the continuum at the far red end of the spectrum. The differences between the model and the data are 
presented as a histogram at the bottom of the plot. 
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Figure 3. The results of minimum fitting to the optical contiimurii of l'KS1549-79 assuming the zero reddening for the nebular 
continuum, restricting the wavelength range of the fits to Xrest < 6OOOA. The plot on the left gives the reduced of the model fits for 
different combinations of YSP age and reddening, while the plot on the right gives the percentage contribution of the YSP component 
in the normalising bin used for the modelling (4720 - 4820A). Reduced values less than unity are deemed acceptable. 



blet adequately using a single Gaussian component that is 
both broad and blueshifted relative to the narrow compo- 
nent detected in the low ionization hnes (see also Tadhunter 
et al. 2001). However, given that results obtained for [OIII] 
in apertures slightly offset from the nucleus show the clear 
presence of a narrow (FWHM < 400 km s"'^) emission line 
component that is also detected in the low ionization lines, 
we have also fitted each component of [OIII] doublet in the 
nucleus using a double Gaussian model, with one of the 
Gaussian components constrained to have the velocity width 
and redshift of the narrow component detected in the low 
ionization lines. The resulting fit — comprising the narrow 
component and an intermediate, blueshifted component — 
is shown in Figure 6. The imperfections of this fit in the line 



cores are to likely refiect the fact that a single Gaussian is 
only an approximation to the true profile of the intermedi- 
ate, blueshifted component in the nucleus. Note that such 
imperfections would not be apparent in fits to the low ioniza- 
tion lines, because of the weakness of the broad components, 
and the dominance of the narrow components in such lines. 

Following on from the double Gaussian fits to [OIII] we 
then attempted to fit all the emission lines in the nuclear 
regions with a double Gaussian model comprising two com- 
ponents: 

- A narrow component (FWHM = 383±15 km s"^ ) with 
redshift 2 = 0.15225, measured from double Gaussian fits to 
the bright low ionisation lines. 
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Figure 4. Example fit to the optical continuum for tiie model including a YSP with age 0.1 Gyr and reddening E{B — V) = 0.5 (zero 
reddening for the nebular continuum). The Call K, G-band and Mgl absorption features are well fitted by this model, but the higher 
order Balmer lines are likely to suffer from residual emission line contamination. 



- An intermediate component (FWHM = 1282±25 km 
s"""^ ) blue shifted by 679±20 km with respect to the 
narrow component, measured from double Gaussian fits to 
the [OIII] lines. 

This model provides an adequate fit to all the lines in the 
nuclear emission spectrum, with the following notable ex- 
ceptions: 

- The Ha+[NII] blend requires an additional broad Ha 
emission line component to fit the wings of the bicud. Wo 
find that we can obtain an adequate fit to the blend if we 
set the width and redshift of the broadest Ha component 
to be the same as those deduced from single Gaussian fits 
to the broad Paa line detected by Bellamy et al. (2003: 
z = 0.15238, FWHM = 1940 km s"^). The flux measured 
for the broadest Ha component is 4.2 x 10~^^ erg cm~^ s~^. 
The detection of this component provides further evidence 
that the reddened quasar nucleus is detected at the longer 
wavelength end of the optical window. 

- The faint high ionization lines [NeV]A3425 and 
[FeVH]A6087 could be fitted adequately by the broad com- 
ponent alone, whereas the faint low ionization [NIJA5200 
could be fitted by the narrow component alone. Presumably, 
in each of these cases the other component that would oth- 
erwise be detected in a double Gaussian fit is not detected 
because of the low S/N and uncertainties in the continuum 
subtraction. 

The velocity shifts and widths of the various kinematic 
components detected in PKS 1549-79 arc compared in Fig- 
ure 5 and Table 2. Not suprisingly, the intermediate and 
narrow components in the double Gaussian fits have similar 
redshifts and line widths to the high and low ionization lines 



Component 



Redshift 



Line width(FWHM) 
(km s-i) 



Nuclear narrow 
Nuclear intermediate 

Paa core 
Paa broadest 

HI 21cm narrow 



0.15225±0.00003 
0.14959±0.00003 



383±15 
1282±25 



0.15238±0.00006 1250±70 
0.1529±0.0002 3525±200 

0.15195±0.00006 80±30 



Table 2. Comparison of the redshifts and line widths of the 
various kinematic components detected in the nuclear regions of 
PKS 1549- 79. Results for nuclear narrow and intermediate com- 
ponents are based on double Gaussian fits to the optical emission 
lines detected in the VLT spectrum, while the Paa results are 
based on double Gaussian fits to the Paa feature detected in the 
new AAT K-band spectrum. 



measured in single Gaussian fits to the lower resolution data 
of Tadhunter et al. (2001). It is notable that both the width 
and velocity shift of the nuclear narrow component are con- 
sistent with those deduced for the extended Ha emission on 
both sides of the nucleus; the redshift of this component is 
also close to that measured for both the core of the broad 
Paa line and the narrow HI absorption (see below). There- 
fore we identify the nuclear narrow component with the rest 
frame of the host galaxy, and the intermediate component 
detected most clearly in the [OIII] lines with a highly dis- 
turbed component that is blueshifted relative to this rest 
frame. 
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Figure 5. Radial velocity profiles (top) and the spatial variations of emission line flux (bottom) for the NTT observations of PKS 
1549-79. The left hand images are for PA-5 and the right hand images are for PA25. Small open and filled triangles represent the narrow 
and intermediate components of Ha respectively. Overplotted is the radial velocity of the deep HI 21cm absorption (large filled triangle 
at -30 km s^-'^ from Morganti ct al. 2001), the shifts of the [O II] and [O III] lines from Tadhuntcr ct al. (2001) (large filled diamonds at 
-5 km s~^ and -550 km s^^ respectively), the 2 components of the [O III] model in the nucleus (large open and filled squares) and the 
broad components of Pa« listed in Table 2 (open 8-pointed stars). For PA-5, positive positions are to the N and negative positions are 
to the S. For PA25, positive positions are to the SSW and negative positions are to the NNE. 



3.2.3 Ionization and physical conditions 

Given the success of the two component emission line profile 
described in the last section, it is possible to investigate the 
ionization, reddening and physical conditions for the nar- 
row and intermediate emission line components detected in 
the nuclear aperture separately. The emission line ratios for 
both of the components derived from the double Gaussian 
profile fits arc shown in the second column of Tabic 3. In 
addition to the nuclear aperture. Table 3 also shows emis- 
sion line ratios determined for a 1.3 x 1.3 arcsecond aperture 
centred 1.6 arcseconds to the east of the nucleus along PA75 



(API: VLT spectrum). In the case of this extended aperture, 
which is affected by seeing disk spillover of the nuclear emis- 
sion line flux, it proved necessary to use a double Gaussian 
fit to model the lines accurately - one component repre- 
senting the nuclear intermediate component and the other 
the extended narrow component. However, despite the con- 
tamination by the nuclear flux in this aperture, the narrow 
component has a much better contrast relative to the inter- 
mediate component than in the nuclear aperture; in Table 
3 we only Ust the narrow line ratios for this aperture. 

It is immediately clear from Table 3 that the 
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Figure 6. Double Gaussian model for the [O III]AA4959,5007 
doublet in the nucleus, forcing the redshift of the narrow compo- 
nent to be consistent with the redshift of the narrow component 
in the low ionisation lines. The composite model is overplotted 
(bold line) on the extracted spectrum (faint line). All of the 4 
components (2 for each line) are also plotted: the dotted lines 
trace the narrow components and the dashed lines trace the broad 
components. 



blueshiftcd, intermediate component detected in the nucleus 
is characterised by a high ioniatiou state, with relatively 
strong [NeV], [Nelll], Hell, [OIII] and [FeVII] lines, and rel- 
atively weak [Oil], [OI] and [SII] lines. 

In contrast, the ionization state of the narrow compo- 
nent detected in a both the nuclear and extended apertures 
is relatively low, with [Oil], [01], [Nil] and [SII] all strong 
relative to the higher ionization [NeV] , [Nelll] and Hell and 
[OIII] lines. 

On the basis of the Ha/H/3 ratio alone there is no strong 
evidence for a high degree of reddening affecting the narrow 
line region in this object, in contrast to the situation in 
PKS1345-I-12 - an object that also shows broad, blueshifted 
Balmer lines (Holt et al. 2003). However, it is difficult to 
be entirely confident of the accuracy of the Ha flux be- 
cause of the complexity of the 7 Gaussian model required to 
fit the Ha-|-[NII] blend. Moreover, the higher order Balmer 
lines cannot be used to check the reddening estimated from 
Ha/H/3 because of their weakness relative to the continuum 
and underlying Balmer absorption line features. Given these 
uncertainties we have made no attempt to correct the emis- 
sion line ratios for reddening. 

Similarly, it has not been possible to estimate ac- 
curate electron temperatures for the narrow line gas us- 
ing the [OIII](5007+4959)/4363 lino ratio, because the 
faint [OIII]A4363 line is highly sensitive to the accu- 
racy of the continuum subtraction, particularly for the 
broader, blueshiftcd component that is also blended with 
H7. Nonetheless it has been possible to estimate the elec- 
tron density from the [SII] (6731/6717) ratio for the narrow 
emission line component detected in the nuclear aperture 
and the extended aperture along PA75. Using a model for 
the [SII]AA6717,6731 blend that includes both an interme- 
diate blueshifted component and a narrow component for 
each line in the doublet, and assuming an electron temper- 
ature of lOjOOOK for the region emitting the [SII] lines, we 



estimate an electron density 430 ± 50 cm^"' for the narrow 
component in the nuclear aperture. Unfortunately, because 
of the relative weakness of the intermediate blueshifted com- 
ponent to the lines, the electron density in the gas emitting 
the blueshifted component is not well-constrained by these 
observations. We find that, although 4 Gaussian models with 
the doublet ratios left unconstrained tend to the high den- 
sity limit for the intermediate component, models with the 
ratio of the intermediate component fixed at the low density 
limit also provide adecjuatc fits to the blend. On the other 
hand, the ratios determined for the narrow components are 
not sensitive to the details of the model used for the inter- 
mediate components. 

In order to investigate the ionization mechanism(s) for 
the various kinematic components, the emission line ratios 
are compared with various ionization models and the results 
for other radio galaxies in the four diagnostic diagrams pre- 
sented in Figure 7. Four sets of models arc plotted: single 
slab optically thick power law photoionization models calcu- 
lated using the MAPPINGS code; mixed medium photoion- 
ization models taken from Binette et al. (1996); pure shock 
cooling zone models and shock plus precursor models from 
Dopita & Sutherland (1995, 1996). 

These diagnostic diagrams further emphasise the high 
ionization character of the gas emitting the intermediate, 
blueshifted component. In all diagnostic diagrams except 
that involving [Nil], the line ratios for this component can 
only be fitted with power-law AGN photoionization models 
that have a high ionization parameter (0.05 < U < 0.1), and 
the line ratios are most consistent with models that have a 
power-law spectral index of -1.5. Although the particular 
mixed medium photoionization models plotted in the dia- 
grams do not provide a good fit to the line ratios of the 
intermediate component, it is possible that, by changing the 
parameters of such models (e.g. ionization parameter, col- 
umn depth of matter bounded component), an acceptable fit 
could be obtained. It is also notable that none of the shock 
or shock-hprecursor can provide an adequate fit to the line 
ratios. 

The intermediate component has line ratios that place 
it beyond the high ionization end of the locus of narrow 
line radio galaxies, overlapping with the range of ionization 
parameter measured for the narrow, blueshifted absorption 
line systems detected against the UV and X-ray continua 
of Seyfert 1 galaxies and quasars (Crenshaw et al. 2003). 
Indeed we would expect such absorption line systems to be 
detected in PKS1549-79 for lines of sight to the AGN that in- 
tercept the high ionization gas. The only line ratio that does 
not fit in with this general trend is [Nil] /Ha, which is un- 
usually large considering the high ionization state revealed 
by the other lines. This may reflect the fact that the system- 
atic uncertainties for this ratio are larger than for the others 
because of the complexity of the Ha-|-[NII] blend. Moreover, 
the models plotted in Figure 7 assume Solar elemental abun- 
dance ratios, whereas the abundance of nitrogen relative to 
other elements such as oxygen may be significantly larger 
than the Solar value - the favoured explanation for the un- 
usually strong [Nil] in, for example, Cygnus A (Tadhunter, 
Mctz & Robinson 1994). 

The diagnostic diagrams also show that the line ratios 
of the narrow emission line components (all apertures) fall 
at the lower ionization end of the ratios measured for the 
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Figure 7. Emission line diagnostic diagrams for PKS 1549-79. Upper left: [OIII]/H/3 vs. [OI]/[OIII]. Upper right: [OIII]/H/3 vs. [NII]/Ha. 
Bottom left: [OII]/[OIII] vs. [NeV]/[NeIII] . Bottom right: [OIII]/[OII] vs. [OI]/[OIII]. In these diagrams, the pure shock models of Dopita 
&; Sutherland (1996) for a range of magnetic parameter {B/^/n = 0, 1, 2 G cm'^/^) with shock velocities 150 < Vnf^^ck 5; 500 km s~^ 
are represented by dotted lines, while the shock plus precursor models are represented by dashed lines. Solid lines represent optically 
thick power law [Fi, oc u"") photoionization models calculated using the MAPPINGS code for three power-law spectral indices {a = -1.0, 
-1.5, -2.0) for a sequence in the ionization parameter (2.5 X 10~* ^ U < 10~^) increasing from right to left. Most of the photoionization 
models plotted in these figures were calculated for low electron densities {N = 10^ cm~'^), however the dot dashed lines shows a sequence 
of models calculated at a single ionization parameter (C/ = 5 X 10~^) for a range of densitities (10^ < AT < 10^ cm-3). The bold line 
corresponds to mixed-medium models from Binette et al. (1996) that include both matter-bounded and ionization-bounded clouds for 
a range of ratios of the solid angles covered by the matter bounded and ionization bounded clouds (10~^ < Aj^.j/j < 10) increasing 
from right to left. The nuclear intermediate component is represented by the large solid square, while the narrow components for the 
nuclear and extended apertures are represented by large solid and open triangles respectively. For comparison, the small filled triangles 
represent data for the nuclear narrow line regions of radio galaxies (data from: Wills et al. 2002, Storchi-Bergmann et al. 1996, Tadhunter 
1987, Grandi & Osterbrock 1978, Costero & Osterbroack 1977, Cohen & Osterbrock 1977, Cohen & Osterbrock 1981, Dickson, 1997, 
Solorzano-Inarrea 2001, Villar-Martin et al. 1998, Clark 1996, Robinson et al. 2000), small filled squares represent extended emission 
line regions (data from: Tadhunter et al. 1994, Storchi-Bergmann et al. 1996, Morganti et al. 1991, Robinson et al. 2000), and open 
squares represent EELRs with strong jet-cloud interactions (data from: Villar-Martin et al. 1998, 1999, Clark 1996, Clarke et al. 1998, 
Solarzano-Innarea 2001). On addition, we present the line ratios for various kinematic components detected in the spatiallty resolved 
NLR of Cygnus A by Taylor, Tadhunter Sz Robinson (2003), with the nuclear broad and narrow components plotted as small open 
triangles, and the extended narrow components plotted as open diamonds. 
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Figure 8. Ncar-IR (K-band) spectrum of PKS1549-79 taken us- 
ing tlie IRIS2 speetrograpli on tlie AAT. Tiie main spectral fea- 
tures are marked. At near-IR wavelengths the spectrum is domi- 
nated by a quasar component. 

narrow line regions of radio galaxies in general, consistent 
with photoionization at moderate ionization parameter, or 
shock-|-precursor ionization models. Despite the evidence for 
young stellar populations in this source, and the possibility 
of on-going star formation, the strengths of the [OI] and 
[Nil] emission lines relative to Ha preclude HII region-like 
photoionization by normal OB stars. 

3.3 Near-IR spectroscopy 

Our new deep K-band spectrum of the nuclear regions taken 
with the ATT/IRIS2 is shown in Figure 8. This represents 
a substantial improvement over the previous K-band spec- 
trum published in Bellamy et al. (2003) , which was affected 
by a faulty grism that resulted in a splitting of all spectral 
features. As well as the strong Paa emission line, the spec- 
trum shows a feature near 2.26/im that is likely to represent 
a blend of H2i^=l-0S(3) and [SiVI]A1.963Mni, an H2Z/=1- 
0S(2) line and (less certainly) a broad feature close to the 
expected wavelength of redshifted Br<5. 

The Paa feature can be modelled £is a combination of 
two Gaussians of width 1220±90 km s"^ and 3540±320 km 
s"^ (FWHM) with fluxes of (1.10 ±0.16) x 10"" erg cm'^ 
s"^ and (1.92 ± 0.15) x 10"" erg cm"^ s"^ respectively; 
redshifts for these components are given in Table 2. Similar 
results are obtained from fitting our NTT/SOFI spectrum. 
Although the narrower of these two broad components has 
a width similar to that of the intermediate component de- 
tected in the [OIII] lines, its redshift is close to that of the 
narrow emission line components (see Table 2), therefore it 
cannot be emitted by the same gaseous system that emits 
the intermediate [OIII] line component. Given this, and the 
fact that the broader Paa component has a width that falls 
outside the normal range for the narrow line gas associated 
with AGN, it is likely that the Paa line is dominated by 
genuine broad line region (BLR) emission. For comparison, 
a single Gaussian fit to the Paa gives a linewidth: FWHM 
= 1940±50 km s~\ 

The new spectra also allow us to resolve the uncer- 
tainty in the identification of the feature at ~2.26/im in 
the observed spectrum (see Figure 8). In terms of its red- 
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Table 3. Data for all optical emission lines detected in the nu- 
clear and extended apertures of PKS 1549-79. Columns are; (a) 
emission line; (b) emission line component where n and i are 
the narrow and intermediate components; (c) emission line fluxes 
measured for a 1.3x 1.3 arcseond aperture centred on the nucleus, 
extracted from the PA75 VLT data (fluxes expressed relative to 
H/3); (d) emission line fluxes measured for a 1.3x1.3 arcseond 
extended aperture centred 1.6 axcseconds east of the nucleus, ex- 
tracted from the PA75 VLT data. 

f weighted mean of [O 11] doublet at the low density limit. 
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Figure 9. Fit to the shape of the near-IR H+K spectrum of PKS1549-79 talcen on the NTT, showing the extrapolation to optical 

wavelengths for various assumptions about the power-law shape of the optical-IR continuum. Also shown are the VLT optical spectrum 
(upper curve at short wavlengths), and the VLT spectrum following subtraction of the best fitting stellar continuum model (lower curve 
at short wavelengths). The upper and lower dashed lines represent models with the minimum and maximum reddening respectively, 
while the solid line represents the average reddening case (see Table 4). The arrow in the lower left corner of the plot represents the 
upper limit on the AGN continuum flux at 5900A derived from HST imaging observations. 



shift this feature could either be identified with blueshifted 
[SiVI]A1.963^im or with H2!^=1-0S(3) at the redshift of the 
nuclear narrow component. However, a single Gaussian fit 
to this feature gives a width of 970±80 km s~^, which is 
significantly broader than the nuclear narrow component, 
but significantly narrower than the intermediate compo- 
nent, suggesting that the feature is composite. Therefore we 
have attempted a double Gaussian fit to the feature, with 
one Gaussian component representing H2t'=l-0S(3) with 
the same width and redshift as the nuclear narrow compo- 
nent, and the other representing [SiVI] A1.963/im with the 
same width and redshift as the nuclear intermediate com- 
ponent. This double Gaussian model provides a good fit to 
the data with fluxes of (4.5 ± 1.5) x 10~^® erg cm~^ and 
(3.1±0.3) X 10"^^ erg cm"^ for the narrow H2i/=1-0S(3) and 
intermediate [SiVI] A1.963/irn components respectively. The 
detection of H2f = 1 — S(3) is consistent with the presence 
of H2i'=l-0S(2) in the spectrum at longer wavelengths; the 
estimated luminosity for the H2i^=l-0S(3) line (2.3x lO'*'' erg 
s^^) is less than estimated by Bellamy et al. (2003) who as- 
sumed that all the ^^2.26/1111 feature is H2!^ = 1 — S(3). It is 
notable that the detection of the blueshifted [SiVI]A1.963/xm 
line is consistent with the high ionization character of the 



intermediate emission line component deduced from the op- 
tical observations (see section 3.2). 

Along with information about the emission line systems, 
the new spectra provide key information about the quasar 
continuum emission. In particular, the new NTT spectrum 
has a broader spectral coverage than the AAT/IRIS spectra 
and allows more stringent limits to be placed on the red- 
dening to the AGN. To estimate the AGN reddening we fol- 
low the same tcchuicjuo as Bellamy ct al. (2003): fitting the 
near-IR continuum spectrum assuming a range of intrinsic 
power-law optical-infrared spectral shapes for quasars taken 
from Simpson & Rawlings (2000). The results are presented 
in Table 4 and Figure 9. 

On the basis of this modelling alone it is clear that 

the V-band absorption and absolute magnitude of the AGN 
rrmst fall in the range 4.9 < A^iintrinsic) < 13.2 and 
—27.56 < My < —22.64 respectively. However, the recent 
HST/ACS imaging of PKS1549-79 (Batcheldor et al. 2006 
in preparation) allows us to place a stringent upper limit on 
thefluxof the quasar core at 5900A(rcst) of Fa < 1.0x10"^'^ 
erg cm^"^ A^^. Using this limit it is clear from Figure 9 that 
the intrinsic absorption and luminosity of the quasar must 
lie at the upper end of the range predicted on the basis of of 
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Minimum 


Average 


Maximum 


« (intrinsic) 


1.62 


0.90 


-0.67 


E(B-V) 


1.52 


2.31 


4.04 




4.94 


7.52 


13.19 


M„ 


-22.64 


-24.18 


-27.56 




-26.97 


-27.43 


-28.44 



Table 4. Reddening and absolute magnitude estimates for the 
quasar nucleus in PKS1549-79. The second, third and fourth 
columns give the minimum, average and maximum reddening cor- 
responding to the range of optical/UV power-law slopes deduced 

by Simpson & Rawlings (2000). The second row gives the values 
of the assumed spectral index for the intrinsic optical-IR spectral 
shape (convention: Fi, oc u^"). The final four rows give estimates 
for the E(B-V) reddening, A„ extinction, and the estimated ab- 
solute magnitude of the quasar nucleus in both the V and the K 
bands. In modelling the continuum shape we have assumed the 
extinction law of Rieke & Lebofsky (1985). 

the ncar-IR continuum modelling; wc find that the optical 
flux limit implies > 6.4, My < -23.5 and Mr < -27.3. 
This places the AGN in PKS 1549-79 firmly in the luminos- 
ity range of quasars {My < —23.0). Note also that such 
high values of nuclear extinction and luminosity are consis- 
tent with the residual spectrum obtained by subtracting the 
best fitting stellar continuum model from our VLT spectrum 
(see Figure 4). 

3.4 Radio imaging and spectroscopy 

3.4-1 VLBI continuum imaging 

Figure 10 shows the VLBI images made with SHEVE at 2.3 
and 8.4GHz. The radio source is characterised by a compact 
core to the western extreme of the source, which shows a 
jet-like extension in PA64, and a highly coUimated jet struc- 
ture to the cast aligned along PA120. The overall extent of 
the source in the 2.3GHz image is ISOmas (430 pc), and 
there is no strong evidence for emission from more extended 
structures; >95% of the total (single dish) flux is contained 
within the structures visible in the VLBI images. Relative 
to the western core the jet appears significantly fainter in 
the 8.4GHz image (core/jet ratio ~10:1) than it does in the 
2.3GHz image (core/jet ratio ~2:1), but the western core has 
a similar brightness in both images. The non-simultaneity of 
the 2.3GHz and 8.4GHz VLBI data, and the resolution dif- 
ferences between the two frequencies, preclude accurate de- 
termination of the power-law spectral indices of the jet and 
core components. However, the data are consistent with a 
steep spectrum for the jet (q ^ -1.2; oc f^") and a 
much flatter spectrum (a ~ 0) for the core. Note that the 
total flux recovered in our VLBI observations at 8.4 GHz 
(3.6 Jy — data taken in 1991) is significantly larger than 
the 8.4 GHz flux measured more recently by Drake et al. 
(2004b: 2.4 Jy — data taken between 1996 and 1998). This 
suggests that the core is variable at high frequencies. 

Overall, the radio structure is similar to the core-jet 
structures detected in some other radio-loud quasars such 



as 3C273. Such structures are generally interpreted in terms 
radio emission from a two sided relativistic jet structure, 
with the axis of jets pointing close to our line of sight. In 
this case the radio emission from the jet on the near side is 
enhanced because the radio emission is beamed towards the 
observer, whereas the emission from the jet on the far side 
is dimmed because the radiation is beamed away from the 
observer. The significant curvature in the jet is also consis- 
tent with the beaming hypothesis, since a small degree of 
intrinsic jet curvature will be exaggerated if the jet is point- 
ing close to the line of sight. If this beaming/orientation 
explanation for the one-sided jet is correct, then the total 
radio emission from the source will be boosted relative to 
the situation in which the jet is close to the plane of the 
sky. Therefore, although the total radio power of PKS1549- 
79 uncorrected for beaming {P2.3GHZ = 5.4 x 10^'' W Hz"^) 
would lead to its classification as a powerful radio galaxy 
that would normally appear as an FRII radio source if ori- 
entated close to the plane of the sky, the true (unbeamcd) 
radio power of the source could be much lower. If we assume 
that the radio flux not recovered in the VLBI maps (i.e. 5% 
of the total) represents an unbeamed radio component, this 
leads to a lower limit on the total unbeamed radio power of 
P2.3GHZ = 2.7 X 10^^ W Hz"^ - just above the break in the 
radio luminosity function (Auriemma et al. 1977) and in the 
transition region between FRI and FRII sources (Fanaroff 
& Riley 1974). 

Assuming that the one sided nature of the jet is a conse- 
quence of beaming effects, we can use the lower limit on the 
jet to counter-jet ratio {Rjc) derived from the radio maps to 
estimate the inclination of the jet to the line of sight. The 
21cm LBA continuum map (see Figure 11) is the most sen- 
sitive for the detected of extended jet features close to the 
nucleus. The peak flux in the jet detected in this map is 2.3 
Jy beam~^, whereas the 3cr detection limit for the counter 
jet is 48 mjy beam"^. Therefore a lower limit on the jet to 
counter-jet ratio is Rjc > 48. The jet to counter-jet ratio 
is related to the bulk Lorentz factor in the jet (7) and the 
inclination of the jet to the line of sight (i) by: 



/ 2 \ 1/2 / pl/{2-a) , \ 



where a is the spectral index of the jet continuum emission. 
For a = —1 and Rjc > 48 we estimate an upper limit on 
the inclination of i < 55 degrees, assuming that the Lorentz 
factor for the bulk motion of the jet falls within the typical 
range estimated for extragalactic radio jets (5 < 7 < 10). 

An alternative possibility is that the one-sided jet is a 
consequence of a stronger interaction between the radio jet 
and the interstellar medium on the east side of the nucleus 
that leads to local re-acceleration of the electrons and en- 
hanced radio emission; a strong jet-cloud interaction would 
also explain the curvature in the jet. However, this scenario 
would not account for the relative strength of the flat spec- 
trum core and the fact that the asymmetric character of jet 
extends into the near-jet close to the core, as is clear from 
the higher resolution 2.3GHz map (Figure 10 middle). 



© 0000 RAS, MNRAS 000, 000-000 



The compact radio source PKS1549-79 15 




l k» ; 

0-1 0-05 -0-05 -0.1 



RA offset (arcsec) 




0.1 0.05 -0.05 -0.1 



RA offset (arcsec) 

Figure 10. High resolution VLBI radio maps of PKS1549-79. 
The upper map is for a frequency of 2.3GHz excluding the Har- 
tebeesthoek antenna, the middle map is also for a frequncy of 
2.3 GHz but including the Hartebeesthoek antenna, while the 
lower map is for an observation frequency of 8.4 GHz (contours: 
-10,10,20,40,80,160,320,640,1280 mJy/Beam). 

3.4-2 HI 21cm kinematics 

Given the core-jet character of the radio continuum struc- 
ture, which suggested that the jet is pointing close to our 
line of sight, it was a surprise to detect HI absorption in 
our ATCA observations (Morganti et al. 2001). On the ba- 
sis of the simplest versions of the unified schemes one would 
normally expect a relatively unobscurod view of the quasar 
nucleus from the direction of the radio jet (Barthel 1989). 
The absorption in the spatially integrated radio emission has 
a FWHM~80 km s~^ and the peak optical depth is about 
2%. This leads to a column density of (4.0 ±0.5) x lO^^Tsp^n 
cm~'^, where Tspin is the HI spin temperature. A possible 
interpretation described in Tadhunter et al. (2001) is that 
the HI absorption is associated with a dusty cocoon around 
the nucleus that also emits the narrow emission line com- 
ponent, and is responsible for the obscuration of the quasar 
nucleus. 

The new LBA data provide information about the spa- 
tial structure of the HI absorption. Figure 11 shows the pro- 



files of the HI absorption measured at the locations of both 
the core and the jet structures. Consistent with the idea of 
a foreground absorbing screen proposed by Morganti et al. 
(2001) and Tadhunter et al. (2001), the HI absorption was 
detected against both of the major radio components with 
similar redshift, linewidth and optical depth. Therefore the 
diameter of the HI absorbing region is > 120 milliarcseconds 
(> 286 pc). Note that the narrow HI absorption is signifi- 
cantly blueshifted (by ~75 km s"'^) relative to the narrow 
optical omission lines (sec Figure 4 and Table 2), but this 
blueshift is an order of magnitude smaller than that mea- 
sured in the [OIII] lines, and might be attributed to gravi- 
tational motions if the absorbing cloud has not yet settled 
into the circumnuclcar disk. 

In Figure 12 the HI absorption profile from the ATCA 
observations (from Morganti et al. 2001) is superimposed 
on the integrated absorption profile obtained from the new 
LBA observations. This shows that most of the HI absorp- 
tion is recovered in the high spatial resolution LBA obser- 
vations, and that the HI kinematics in the LBA and ATCA 
are consistent. 



4 DISCUSSION 

4.1 Evolutionary status and the triggering of the 
activity 

The ready detection of multiple tidal tails at a relatively 
high surface brightness level is entirely consistent with the 
idea that the activity in PKS1549-79 has been triggered in 
a major merger between galaxies, at least one of which was 
gas rich and had a cold disk component prior to the merger. 
Indeed, the tail structure observed in PKS 1549-79, with two 
near-parallel tails pointing to the south and a jet-like struc- 
ture pointing to the north bears marked similarity with the 
structures produced in numerical simulations by Barnes & 
Hernquist (1996) of the merger between two disk galaxies, 
observed from the orbital plane of the merging galaxies. In 
such a merger, it will take of order 500 Myr for the nuclei 
of the two galaxies to finally merge together. Therefore it is 
likely that that we are observing the system a considerable 
period after the start of the encounter. 

The young stellar populations (YSP) provide useful 
complementary information about the merger. The presence 
of a reddened YSP in the nuclear regions of this object is 
consistent with the idea that the merger also triggered a 
starburst. Given that the YSP has an estimated age of 50 
- 250 Myr at the observation epoch, it is likely that the 
starburst was triggered a considerable period after the start 
of the merger sequence. This in turn suggests that one of 
the galajcies had a significant bulge that hindered the for- 
mation of a bar that would otherwise lead to the major 
merger-induced star formation occuring early in the merger 
(Mihos & Hearnquist 1996). Notliwithstanding the uncer- 
tainties about the heating mechanism for the dust emitting 
the far-IR radiation, the substantial far-IR excess in this 
source suggests the presence of substantial ongoing star for- 
mation in the nuclear regions that is simultaneous with the 
currently observed quasar and jet activity. 

It is also interesting to consider at what stage the rar 
dio source was triggered in the merger sequence. Tadhunter 
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Figure 11. Middle: 1.64 GHz LBA continuum image of PKS1549-79 obtained from the line-free channels. Right: integrated profile of 
the HI absorption in the western core source. Le/i;: integrated profile of the HI absorption for the eastern jet source. 



et al. (2001) have argued that, since the extreme forbidden 
emission line widths measured in PKS1549-79 are only de- 
tected in radio galaxies with young, compact radio sources, it 
is likely that we are observing the radio source at a relatively 
early stage after the triggering of the activity (< 10''yr). If 
this is correct, it is plausible that there was a significant 
delay (>40 Myr) between the starburst associated with the 
YSP detected at optical wavelengths and the triggering of 
the radio jet activity. However, the evidence for such a time 
delay is weaker in PKS1549-79 than for some other radio 
sources in which the time lag between the the major episode 
of merger-induced star formation and the triggering of the 
currently observed radio source is estimated as 0.3 ~ 2.0 Gyr 
(Tadhunter et al. 2005, Emonts et al. 2005, Johnston et al. 
2005). 



4.2 The nature of the AGN 

Our new data provide key information about the nature of 
the AGN in PKS1549-79. Most notably, the relatively nar- 
row Paa detected in the quasar component of PKS1549-79 
(FWHM = 1940±50 km s~^ from a single Gaussian fit) is 
consistent with the classification of this source as a narrow 
line Seyfert 1 galaxy (NLSl: FWHM < 2000 km s'^ Oster- 
brock & Pogge 1985, Osterbrock 1987) or more accurately, 
given its luminosity, a narrow line quasar. Interestingly, the 
extreme blueshift of the [OIII] line relative to the galaxy 
rest frame in this object is also similar to those detected 
in the broad [OIII] emission line components of many NLSl 
galaxies and similar objects (Veron-Cetty & Veron 2001, Za- 
manov et al. 2002, Boroson 2005, Aoki, Kawaguchi& Ohta 
2005). 

NLSl are often interpreted as AGN that are accreting 
at a relatively high Eddington radio {Ltoi/ Ledd ~ 1), based 
on the small virial black hole mass estimates deduced from 
their relatively narrow permitted line widths (Pounds, Done 
& Osborne 1995, Mathur 2000, Mathur & Grupe 2005). It 
has also been suggested that the NLSl fall significantly be- 
low the correlations between black hole mass and host galaxy 



VLBI + ATCA 




Figure 12. Profile of the HI absorption in PKS 1549-79. The 
solid line represents the ATCA data of Morganti et al. (2001), 
while the dotted lines represent the integrated HI profile measured 
from the new LBA data. Most of the absorbed flux is recovered 
in the higher resolution LBA observations. 



bulge properties (Grupe & Mathur 2004). However, it is im- 
portant to add the caveats that (a) the small permitted line 
widths of NLSl may reflect orientation effects (e.g. BLR 
disk viewed face on) rather than low black hole masses; and 
(b) the bulge properties for NLSl in the Grupe & Mathur 
(2004) study were deduced from [OIII] linewidths, which are 
likely to be affected by non-gravitational (radial) motions 
suggested by the large [OIII] blueshifts detected in some 
sources. 

It is interesting to consider whether PKS1549-79 is also 
a high accretion rate object with a large Eddington ra- 
tio. Its black hole mass can be estimated using two inde- 
pendent methods. First, a virial mass estimate can be ob- 
tained using the Paa line width and luminosity. Assuming 
the mimimum extinction for Paa {Ak = 0.74 magnitudes). 
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and a Paa/H/3 ratio set at the case B recombination value 
(Paa/ap = 0.332, see Osterbrock 1989), we deduce an in- 
trinsic quasar H/3 luminosity of L(Hp) = 1.0 x 10''^ erg 
s~^ for our adopted cosmology. Assuming that FWHM(H/3) 
= FWHM(Paa) = 1940 km s"^ wc then use the relation- 
ship between FWHM(H/3), L(H/3) and M^h published by 
Greene & Ho (2005: their eqn 7) to estimate a virial black 
hole mass of Mbh = 5.3 x lO'^ Mq. A second black hole 
mass estimate can be obtained from the correlation between 
R-band bulge luminosity (A^h) and black hole mass pub- 
lished by McClure and Dunlop (2002) along with the pub- 
lished R-band absolute magnitude for PKS1549-79 (Drake 
et al. 2004a: Mr, = -21.9 for our cosmology)*. This leads 
to Mbh = 2.6 X 10^ M0 — almost an order of magnitude 
larger than the virial estimate of black hole mass. However, 
the R-band luminosity estimate can be regarded as an up- 
per limit because the YSP in the bulge and tail features of 
PKS1549-79 along with any AGN contamination are likely 
to artificially boost the R-band luminosity; the only way in 
which the R-band magnitude can substantially underesti- 
mate the black hole mass is if the host galaxy as a whole is 
subject to substantial intrinsic extinction. In contrast, the 
virial estimate is likely to represent a lower limit on the black 
hole mass because of the effects of orientation on the Paa 
line width, especially given that the radio observations sug- 
gest that radio jet is pointing close to the line of sight, with 
the putative BLR disk observed close to face-on. The true 
mass of the black hole in PKS1549-79 is therefore likely to 
fall between the virial and bulge luminosity estimates. 

The other element required to determine the Edding- 
ton ratio is the bolometric luminosity {Lboi) of the quasar. 
Lboi can be estimated using the quasar continuum modelling 
results along with information on generic quasar spectral en- 
ergy distributions (SEDs) presented in Elvis et al. (1994). 
The range of optical SEDs that are consistent with our op- 
tical/IR spectroscopy and imaging have monchromatic V- 
band luminosities in the range 1.2 x 10"'^ < Lv < 4.8 x 10"'^ 
erg s"^ A"^ or 6.6 x lO''* < \Lx < 2.6 x lO"® erg s"^ 
(A = 5500A) for the quasar in PKS1549-79. The V-band 
bolometric correction factor listed in Elvis et al. (1994) 
is: Lboi/i^Lv) = 14.2. Therefore, we estimate that the 
bolometric luminosity of PKS1549-79 falls in the range 
9.0 X 10"*' < Lboi < 3.8 X 10"^ erg s"\ Note that, although 
the radio jet in this source is pointing close to our line of 
sight, there is no evidence that the optical/IR continuum is 
substantially affected by contamination by a beamed syn- 
chrotron source: not only does the spectral energy distri- 
bution of PKS1549-79 show a steep decline between radio 
and sub-mm wavelengths (compare Drake et al. 2004b with 
Beasley et al. 1997), but both the H/? equivalent width and 
the Lboi/Lir ratio fall well within the range measured for 
un-beamed quasars. 

These estimates of Lboi lead to Eddington ratios in the 
range 0.3 < Lboi/L^dd < H for our upper estimate of the 

* Note that the R-band absolute magnitude derived by Drake et 
al. (2004a) is based in the integrated light measured within the 
outermost ellipse fitted to their R-band image (corresponding to a 
semi-major axis of 12.5 arcseonds or 30 kpc). Since this magnitude 
corresponds to the total light, and no disk/bulge separation has 
been attempted, it is likely to overestimate the luminosity of the 
bulge. 



black hole mass {Mbh = 2.6 x 10* Mq) and in the range 
1.3 < Lboi/Ledd < 57 for our lower estimate of the black 
hole mass {Mbh = 5.3 x lO'^ M©). The measurement of such 
high Eddington ratios reinforces the link between PKS1549- 
79 and NLSl galaxies. 

Although we cannot rule out the idea that PKS1549-79 
is a significant outlier relative to the statistical correlations 
we have used to derive both black hole mass and bolomet- 
ric luminosity estimates, our results are at least consistent 
with accretion at a high Eddington ratio in this object. On 
this evidence, accretion at high Eddington ratio does not 
prevent the formation of powerful, relativistic jets. For com- 
parison, the Eddington ratio measured for Cygrnis A the 
only quasar with a direct, dynamical black hole mass esti- 
mate - is (1.6 - 6.0) X 10"^ (Tadhunter et al. 2003). Simi- 
larly, Dunlop et al. (2003) derive relatively small Eddington 
ratios (typically Lboi/Ledd < 0.1) based on R-band black 
hole mass estimates of several radio-loud and radio-quiet 
quasars. On the other hand, Punsly & Tingay (2005) have 
recently highlighted the high redshift quasar PKS0743-67 as 
another object that is accreting at high Eddington ratio but 
has powerful radio jets. 

4.3 Characterising the warm outflow 

Recent theoretical studies of galaxy evolution have invoked 
quasar-induced outflows to explain both the correlations be- 
tween the masses of black holes and the galaxy bulge prop- 
erties (e.g. Silk & Rees 1998, Fabian 1999, di Matteo et al. 
2005) , and the form of the galaxy luminosity function at the 
high luminosity end (e.g. Benson et al. 2003). Such feedback 
models require that a relatively large fraction of the avail- 
able accretion power of the quasars is thermally coupled to 
the circumnuclear ISM (>10%: Fabian 1989, di Matteo et 
al. 2001). 

PKS1549-79 appears to represent a situation close to 
that modelled in some of the recent simulations (e.g. di Mat- 
teo et al. 2005), in which the black hole is growing rapidly 
via merger-induced accretion following a merger between gas 
rich galaxies. We have also found - in the blueshifted [OIII] 
lines - clear evidence for an outflow in the central regions of 
this source. Therefore it is important to investigate whether 
the properties of the waxm gas outflow in PKS1549-79 are 
consistent with the assumptions of the quasar feedback mod- 
els. 

Wc start by determining the mass outflow rate in the 
intermediate (blue shifted) component. For a steady state 
outflow the mass outflow rate is given by: 

M = NeVoutArUp (2) 

where is the electron density, e the volume filling factor, 
Vout the outfiow velocity, A the area of the outfiow and nip 
the mass of the proton. The filling factor is related to the 
H/9 luminosity {L{H/3)), density and emitting volume V as 
follows: 

e=_M_ (3) 

where a'j^Q is the effective H,5 recombination coefficient 
(see Osterbrock 1989), vh/} the frequency of H/3, and h the 
Planck constant. It follows that the mass outflow rate can 
be written as: 
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and for a spherical outflow this reduces to: 
^^SL^^rn^ (5) 

where r is the radius of the spherical volume. On the ba- 
sis of our spectroscopy we estimate an H/3 luminosity of 
L{H0) — 2.8 X 10*" erg s"^ for the intermediate, blueshifted 
component and a rest frame outflow velocity of Vout = 680 
km s"\ If we assume case B recombination theory for an 
electron temperature T = 10, 000 K and that the spherical 
outflow has the same radius as the radio source (180 milli- 
arcseconds or r = 430pc for our cosmology) wo find that the 
mass outflow rate in units of solar masses per year is given 
by: 

M = 9.l(H!^) M^yr-\ (6) 

The largest uncertainty in this calculation is the electron 
density of the outflowing gas (A^). Because it has not been 
possible to measure the electron density for the outflow- 
ing gas directly in the case of PKS 1549-79, we assume that 
it falls in the range measured for the NLR of other radio 
galaxies. The typical electron densities measured for kpc- 
scale ionized gas in extended radio sources such as Cygnus 
A (Taylor, Tadlmnter & Robinson 2003) and 3C321 (Robin- 
son et al. 2000) are N ~ 10^ cm~^. However, the densities 
may be larger in the emission line regions associated with 
compact radio sources, and in the case of PKS1345+12 
a similar object that has broad, blueshifted [OIII] lines - a 
relatively high electron density (A'^ > 5 x 10^ cm~*) has been 
deduced for the blueshifted emission line components (Holt 
et al. 2003). Therefore it is reeisonable to assume densities 
in the range 10^ <N < 10" cm , leading to mass outflow 
rates 0.12 < M < 12 M0 yr~^. It is notable that this covers 
same range as the mass outflow rates deduced for the NLR 
of Seyfert galaxies (Veilleux, Cecil & Bland-Hawthorn 2005) 
and the high ionization absorption line systems detected at 
UV and X-ray wavelengths in some AGN (Crenshaw et al. 
2003) , but is somewhat less than deduced for the neutral out- 
flows detected in extreme starburst sources such as ULIRGs 
(Rupke et al. 2005a, b). For the same range of densities we 
estimate that the total mass in the warm gas outflow falls in 
the range 1.9 X 10* < Mtotai < 1.9 x 10® M© and the filling 
factor falls in the range 2.3 x 10"^ < e < 2.3 x 10"^. 

EYom the point of view of gauging the likely impact of 
the warm gas outfiow on the circumnuclear gas in the bulge 
of the host galaxy of PKS1549-79 it is also important to es- 
timate the kinetic power of the outflow, including both the 
radial and turbulent components in the gas motion. Assum- 
ing that the relatively large linewdith of the outflowing gas 
reflects a turbulent motion that is present at all locations in 
the outflow region, the kinetic power is: 

E = 6.34 X lQ^^^-(vl^t -h (FVKHM)V5.55) erg s"^ (7) 

where FWHM is the full width at half maximum of the 
blueshifted [GUI] line in km s~^, and M is the mass outflow 
rate expressed in solar masses per year. For Vout = 680 km 
s"\ FWHM = 1282 km s"^ and the mass outflow rates 



determined above, we find that the kinetic power for the 
outflowing gas falls in the range 2.2 x 10*° < < 2.2 x 
10*^ erg s~^. Taking the Eddington luminosity derived from 
the larger black hole mass estimate (see section 4.2 above) 
as an estimate of the total power available from accretion, 
the kinetic power in the warm gas component represents 
only a small fraction (7 x 10"'' < E/Ledd <7x 10~^) of the 
available accretion power, in contrast to the quasar feedback 
models which require a much larger fraction of the accretion 
power of the black hole to be thermally coupled to the gas 
surrounding the black hole. It is also unlikely that the warm 
gas outflow by itself is capable of removing all the warm/cool 
gas from the central regions of the host galaxy. Assuming 
that PKS1549-79 has a total mass Mtotai ~ lO" M©, and a 
gas mass Mgas ~ 10^° M© contained within a radius of 5 kpc 

- conservative assumptions for ULIRGs which are likely to 
have a larger gas mass concentrated within a smaller radius 

- the gravitational binding energy of the gas is Ebind ~ 
GMgasMtotai / Rgas ~ 2 X 10"''** Org. In comparisou, the warm 
gas outflow will deposit only ~ 10^^ - lO^'^ erg into the 
surrounding ISM, assuming that it can persist in its current 
form for the typical 10^ yr lifetime of an extragalactic radio 
source. 

How do we explain the fact that the estimated kine- 
matic power in the warm gas outflow in PKS1549-79 is at 
least three orders of magnitude less than might be expected 
on the basis of the quasar feedback models? Considering first 
the assumptions that have gone into the above estimates, 
the key parameters are the radius of the outflow, the elec- 
tron density and the emission line luminosity. E would be 
underestimated if the radius and/or density of the outflow 
have been overestimated, or the H/3 luminosity has been un- 
derestimated. However, it is unlikely that the radius of the 
outflow region has been overestimated, because recent HST 
imaging observations of PKS1549-79 demonstrate that it has 
a similar scale to the radio source (Batcheldor et al. 2006); 
allowing for projection effects, it is more likely that the ra- 
dius we have used is conservative and the outflow region is 
significantly larger than we have estimated. Similarly, given 
that warm gas has been accelerated to large velocities and 
other compact radio sources show strong evidence for a rel- 
atively high density NLR (Holt et al. 2003, Holt 2005), it 
is unlikely that the electron density in the outflow can be 
substantially less than the measured density in the narrow 
component (430 ± 50 cm"^) and the lower hmit {N > 10^ 
cm^"') we have assumed. 

Given the fact that the quasar nucleus in PKS1549-79 
is itself substantially obscured at optical wavelengths, it is 
less easy to dismiss the idea that a portion of the outflow 
region is highly extinguished at optical wavelengths so that 
the H/3 luminosity is underestimated. However, it is unlikely 
that the H/3 has been underestimated by more than a factor 
50 because the blueshifted outflow component would then 
make a substantial contribution to the Paa flux, which is 
not observed. 

Most plausibly, the warm gas outflow we have detected 
in [QUI] represents only a small fraction of the total mass 
in the outflow, with the remaining mass locked in hotter or 
cooler phases of the ISM that are difficult to detect at optical 
wavelengths. Recently, Morganti et al. (2003,2005a,b) have 
detected massive neutral outflows in several compact radio 
sources using observations of the 21cm line of neutral hy- 
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drogen. The mass outflow rates and energy fluxes contained 
within such neutral outflows are substantiaUy larger than we 
have estimated for the waxm gas outflow in PKS1549-79. It 
is notable that there is a hint of a blue wing to the HI 21cm 
line detected in PKS1549-79 (sec Figures 11 and 12). If this 
can be confirmed with higher quality observations, it may 
be possible to make a more complete census of the warm 
and cool gas outflows in this object. 

4.4 HI spin temperature 

The exact value of the HI spin temperature {Tspin ) is a ma- 
jor uncertainty in calculating the column density of neutral 
hydrogen using HI 21cm absorption line observations. Al- 
though Tspin is often assumed to have values consistent with 
those estimated in the Galaxy (100 < Tspin < 1000 K: Heiles 
& Kulkarni 1988), higher values may be appropriate in the 
central regions around AGN because of excitation of the up- 
per level of the 21cm hyperfine transition by collisions, Lya 
photons or 21cm continuum photons (see Bahcall & Ekers 
1969 for a discussion). The latter may be particularly impor- 
tant in a case like PKS 1549-79 because the radio emission 
is beamed towards the HI absorbing cloud along our line 
of sight. However, to our knowledge there have so far been 
no direct estimates of the HI 21cm spin temperature in the 
vicinity of luminous AGN. 

PKS1549-79 provides an excellent opportunity to es- 
timate the spin temperature directly because we have an 
independent estimate of the neutral column towards the 
quasar nucleus from the extinction estimates derived from 
our quasar continuum modelling. The SED modelling gives 
2.0 < E{B — V) < 4.0 for the intrinsic line of sight reddening 
towards the quasar nucleus. Assuming the standard Galac- 
tic dust to gas ratio {E{B - V)/Nhi = 1.7 x 10"^^ mag 
cm~^: Bohlin et al. 1978) this translates into neutral HI 
columns in the range 1.2 X 10^2 < Nhi < 2.4 X 10^2 cm 
On the basis of our HI 21cm absorption line measurements, 
the neutral HI column is related to the spin temperature by: 
Nhi = (4.0±0.5)x 10^* Tspin cm"^. Therefore, by combining 
the estimates of the neutral HI column derived from both 
the HI 21cm and the quasar SED measurements we find 
that the electron spin temperature must fall in the range 
3000 < Tspin < 6000 K. These measurements provide clear 
evidence that the HI spin temperature in PKS1549-79 is 
higher than the typical Galactic value. 

Implicit in our estimate of Tgpin is the assumption that 
the column reddening the quasar is the same as that causing 
the HI 21cm absorption. Although the quasar nucleus can be 
considered point-like, while the radio omission is clearly spa- 
tially extended in our VLBI observations, our LB A observa- 
tions show that the HI absorption covers the entire spatially- 
resolved radio source, including the flat spectrum core to the 
western jet of the source. Therefore it is entirely plausible 
that the HI 21cm and quasar reddening columns arc asso- 
ciated with the same cloud. Another major assumption we 
have made is that the dust-to-gas ratio in the absorbing 
cloud has the Galactic value. This is called into question 
by recent observations that provide evidence for dust-to-gas 
ratios in the vicinity of AGN that are a factor ~10 or more 
lower than the Galactic value (Maiolino et al. 2001). How- 
ever, if such low dust-to-gas ratios also hold in PKS1549-79, 
the resulting spin temperature would be even larger than 



we have estimated above, further strengthening the evidence 
for large Tspin in this object. In view of this the lower limit 
we have derived for the spin temperature in PKS1549-79 
{Tspin > 3.0 X 10^ K) is likely to be conservative. 



5 CONCLUSIONS AND FURTHER WORK 

The new observations presented in this paper confirm that 
PKS1549-79 represents a luminous quasar observed at an 
early stage of its evolution. The quasar and jet activity have 
been triggered in a major merger which has also triggered 
substantial star formation. Accreting at close to the Edding- 
ton rate, the black hole in this system is growing rapidly and 
simultaneously driving relativistic jets. The growth phase 
of the black hole is substantially obscured at optical wave- 
lengths in this object, consistent with recent galaxy evolu- 
tion models (Hopkins et al. 2005). However, although the 
optical emission lines provide evidence for a high ionization, 
warm outfiow associated with the AGN, the warm outfiow 
currently observed at optical wavelengths is not in itself suf- 
ficient to remove all the gas from the central regions of the 
galaxy. 

PKS1549-79 is a key object because it links together 
several classes of active galaxies, including radio galaxies, 
quasars, ULIRGs, NLSl and AGN with high ionization ab- 
sorption line systems. Moreover it shows many properties 
in common with the recently discovered population of high 
redshift AGN detected at sub-mm wavelengths (Alexander 
et al. 2005). The common thread is that all these classes 
of object are associated with the triggering and feedback 
effects of major galaxy mergers. Although much attention 
is paid to the links between star formation and activity in 
objects at high redshift, this study of PKS1549-79 serves to 
emphasise that some galaxies arc evolving rapidly in the lo- 
cal universe, and that it is possible to use such objects to 
study the coevolution of black holes and galaxy bulges in 
depth. 

Several further observational studies of PKS1549-79 
would enhance our understanding of the relationship(s) be- 
tween AGN activity and galaxy evolution. They include the 
following. 

- Near-IR measurements of the stellar bulge luminosity 
and velocity dispersion to improve estimates of the black 
hole mass. 

- High quality optical spectroscopy observations of a 
wider range of emission line diagnostics to measure the elec- 
tron density accurately, and thereby improve our determi- 
nation of the properties of the warm gas outflow. 

- Deep X-ray spectroscopy observations to refine esti- 
mates of the quasar bolometric luminosity and absorbing 
HI column. Such observations would also allow the links be- 
tween PKS1549-79, NLSl and sub-mm galaxies to be further 
investigated. 

- Multi-epoch VLBI observations to search for superlu- 
minal motions in the jet, and accurately determine the in- 
clination of the jet to the line of sight. 

- Deeper HI 21cm measurements and high resolution op- 
tical observations of the NalD line to investigate the impor- 
tance of neutral gas outflows. 
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